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I. INTRODUCTION 


This progress report summarizes our continuing research into the quantitative non- 
destructive evaluation of composite materials. In previous reports we have described 
investigations of impact damage in graphite-epoxy composites. In this report we shift 
our emphasis somewhat and describe investigations we have carried out to characterize 
porosity in composite media using ultrasonic waves. These investigations are described 
in Section II. Previous reports have also described our investigations of local approxima- 
tions to the Kramers-Kronig relations. In Section III of this Progress Report we present a 
derivation of these relations and show that they may also be applicable to systems that 
could conceivably exhibit considerable dispersion such as composite laminates contain- 
ing porosity. 


n CHARACTERIZATION OF POROSITY IN COMPOSITE LAMINATES 
Ha. INTRODUCTION 

The detrimental effects of porosity on material strength are well known. The work 
of Rose, Hsu, and Adler provides a means of estimating the volume fraction of pores and 
the average pore radius in isotropic elastic media from the value of frequency at which 
the attenuation coefficient becomes frequency independent and the magnitude of the 
attenuation coefficient at that plateau. 1 Quantitative results for the isotropic case depend 
on numerical factors obtained by Gubematis et al. which are functions of the ratio of the 
transverse to longitudinal sound velocities. 2 Adler, Rose, and, Mobley have tested these 
theories by making measurements of attenuation covering a frequency range that 
extended well into the frequency independent plateau. 3 Recent investigations suggest the 
feasibility of applying these methods to estimate the volume fraction of porosity in com- 
posite laminates. 4 ' 10 The experimental results of these investigators suggest that the 
theoretical results obtained by Rose et al. are qualitatively correct even though some of 
the features of wave propagation in layered, anisotropic media are not explicitly incor- 
porated into the scattering model. 

In this report, we consider an approach to deal with the case of anisotropic media, in 
which the magnitude of the attenuation may preclude making measurements at 
sufficiently high frequencies to reach the plateau region. Approaches to materials char- 
acterization based on the frequency dependence (slope) of attenuation are widely 
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employed in medical ultrasonics. 11 A well-known experimental result for frequencies 
lower than that corresponding to the plateau (i.e., / < f high ) but above the Rayleigh 
region (f >f [ow ) is that the attenuation coefficient exhibits an approximately linear 
increase with frequency over a range from f ^ to fhigh • Typically the useful bandwidth 
of the measurement is substantially smaller than this range and falls at an unknown loca- 
tion between these limits. Fortunately, the expression for the attenuation coefficient con- 
tains numerical factors which depend only weakly on the ratio of f high to f iow , varying 
only by a factor of 3 for the ratio of fhigh 10 flow ranging from 10 to 1 to 10,000 to 1. In 
the case of composites with complex lay-ups for which a detailed theory that describes 
the effects of porosity on attenuation may not be available, empirical knowledge of these 
numerical factors obtained from laminates of known porosity might provide an approach 
for estimating the porosity from ultrasonic measurements of similar composites. 

To evaluate this approach we measured the slope of attenuation as a function of fre- 
quency in a set of 5 glass-fiber/epoxy-matrix test specimens with simulated porosity 
(glass spheres) ranging from 0% to 12% (volume fraction) and a set of 5 uniaxial 
graphite-fiber/epoxy-matrix specimens with simulated porosity (glass spheres) ranging 
from 1% to 8%. Good correlation was obtained between the measured slopes and poros- 
ity in each case, suggesting that semi-quantitative estimates of porosity can be achieved 
without measurements in the plateau region and without a quantitative theory. 

Hb. EXPERIMENTAL METHODS 

1) SAMPLE PREPARATION 

The effects of porosity were simulated using solid glass beads, 75 to 150 microns in 
diameter, in 16 ply uniaxial graphite-fiber/epoxy-matrix composites. These composites 
were fabricated using 5208-T300 prepreg tape. Measured amounts of glass beads were 
introduced between the 12th and 13th layers during the lay-up of a 12 by 16 inch compo- 
site. The beads were dusted onto circular regions 2 inches in diameter at sites on a 
square grid with centers 4 inches apart. The sample was autoclaved and cured in an oven 
using a standard cure protocol. The 12 by 16 inch sample was cut into smaller samples 
(approximately 3.75" by 3.75") so that each contained a single zone of "porosity" with a 
volume fraction of 1%, 2%, 4%, 6%, or 8%. 

2) MEASUREMENT METHODS 

The signal loss was measured in transmission mode with a specimen placed in the 
overlapping focal zones of a matched pair of 25 MHz center frequency, 0.25 inch 
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Figure 1. Block diagram of the transmission mode acquisition system used for measurements 
of attenuation in porous graphite/epoxy composites. 
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diameter, 1 inch focal length transducers. Each sample was scanned on a 21 by 21 grid 
in 1 mm steps and the acquired frequency spectra averaged to reduce the effects of spa- 
tial variations of "porosity" within the samples. 

The measurement system used for data acquisition is illustrated schematically in 
Fig. (1). The transmitting and receiving transducers were oriented so that the insonifying 
beam was perpendicular to the surfaces of the sample and were aligned by viewing the 
received signal on a spectrum analyzer. A Metrotek MP215 wideband pulser was used to 
drive the transmitting transducer. The output of a MR106 wideband receiver was routed 
to a stepless gate and the 1.5 ) isec gated signal was subsequently used as the input to the 
spectrum analyzer. A DEC PDP 1 1/73 running the UNIX operating system was used to 
control the motor driven apparatus on a C-scan tank (in which the samples were placed 
for data acquisition) as well as to acquire the data from the spectrum analyzer for storage 
and subsequent analysis. 

He. DATA ANALYSIS 

The signal loss through the composite laminate was obtained by normalizing the 
averaged frequency spectrum with a calibration spectrum obtained from a water-only- 
path trace, 

Signal Loss = log [calibration spectrum ] - log [sample spectrum ] . (1) 

This method of log spectral subtraction is performed to deconvolve effects arising from 
the electromechanical response of the transducers and front-end electronics from the 
sample spectrum. The normalized data were analyzed by performing a Taylor expansion 
around the center frequency / of the useful bandwidth, 

Signal loss -Kq + K^xIJ -f ) (2) 

where K 0 is an estimate of the average signal loss over the useful bandwidth, and K x is 
the rate of change of the signal loss with respect to frequency. This procedure is illus- 
trated in Fig. (2) where the signal loss of a typical spectrum is plotted as a function of 
frequency along with the appropriate Taylor expansion. The usable bandwidth of this 
system was 3 to 12 MHz, where the upper limit was due to the increasing attenuation 
coefficient as function of frequency exhibited by the composite laminates. 

He. METHOD VALIDATION 

In order to validate our methods, measurements were made on a set of 5 glass- 
fiber/epoxy-matrix composites containing controlled amounts of simulated "porosity" 
(0%, 1%, 3%, 6%, or 12% volume fraction). The glass fibers were approximately 12 pm 
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Method Validation 



Figure 3. Validation of the proposed method. Preliminary experiments were performed on 
glass-fiber/epoxy-matrix composites with controlled amounts of "porosity" (simulated by solid 
glass inclusions). 
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in diameter and four to five centimeters in length. The fibers (p = 2.43 ± 0.09 gm /cm 3 ) 
were layed-up by hand in an epoxy resin matrix (p = 1.10 ± 0.01 gm/cm 3 ). Porosity was 
simulated by the random inclusion of solid lead-glass spheres (p = 2.47 ± 0.04 gm/cm 3 ) 
drawn from a distribution with radii ranging between 37 (im and 75 pm. These test sam- 
ples were fabricated with a fiber volume fraction of approximately 8%. One sample was 
fabricated without glass beads in order to serve as a control. 

The results of this control study are presented in Fig. (3) in which the slope of the 
attenuation is plotted versus volume fraction of "porosity". The correlation coefficient 
obtained by performing a linear regression between slope of attenuation and the volume 
fraction of "porosity" is 0.9, suggesting the potential of the method. 

Hd. THEORY 

Rose has derived an expression for the concentration of hollow spheres embedded 
in an elastic medium, 

Cone. = —j— j dk (3) 

3A 2 7t J 0 k 2 

where a (k) is the excess attenuation due to the scattering of sound waves by the pores. 12 
That is, a(k ) represents the increase in the attenuation produced by the addition of hol- 
low spherical pores to the material above that due to the attenuation from the elastic 
background medium. We consider the case in which the excess attenuation due to the 
presence of porosity can be approximated by 


a(f ) = ( excess attenuation ) = (excess K{)x 


- —f 4 f <floW 

J low 

< f flow <f< /high ( 4 ) 

f high f ^ f high • 


as illustrated in Fig. (4). 

This model implies an / 4 dependent Rayleigh scattering in the low frequency 
region, and a frequency independent geometrical scattering in the high frequency region. 
The approximately linear dependence on frequency in the intermediate region is well 
supported from experimental measurements performed on composite laminates. Because 
of the finite usable bandwidth of the experimental apparatus, both f high an d flow 316 usu- 
ally unknown. Nevertheless, the relationship between the volume concentration of poros- 
ity and excess K x is only weakly dependent upon the ratio of f high to f iow . 
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Model Attenuation Coefficient 
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Figure 4. Illustration of the phenomenological model. 
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Figure 5 
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Figure 5. Plot illustrating the slow variation of the Numerical Factor as a function of the 

ratioof W o fo* 




- 11 - 


Inserting Eq.(4) into Eq.(3) leads to 


4 V 

Cone. = — — x — x ( excess AT,) x< 
3 A 2 ^ 2k 


J low 0 flow J f high j 


or 


4 V 

X X 

i +to 

fhigh 

+ i 

3 A 2 ^ 2n 

3 

f low 



It is convenient to express this result as 

Cone. = ( excess K{)x Velocity x [Numerical Factor] , 


( 6 ) 


(7) 


where 


Numerical Factor = 


4 1 

— X — X 

3 A 2 7t 27C 


i + ,„ 

fhigh 

+ i 

3 

f low 



( 8 ) 


The Numerical Factor is a slowly varying function of the ratio / high 1 flow over several 
decades, as illustrated in Fig. (5). Because f high and are rarely known, this slow 
variation is the basis for the predictive usefulness of the model. 


He. RESULTS and CONCLUSIONS 


The results obtained for the graphite-fiber/epoxy-matrix specimens are displayed in 
Fig. (6). The scatter plot displays the resultant slope of attenuation (from the spatially- 
averaged normalized spectrum) for each of the five samples versus the volume fraction of 
"porosity" for that sample. The slope of attenuation correlates well (r = 0.9) with concen- 
tration of "porosity". An examination of Equation(7) indicates that the value of the 
Numerical Factor for the composite laminates under investigation can be determined 
from the slope of the "Slope of Attenuation" versus the "Volume Fraction of Porosity" 
correlation plot, 


Numerical Factor = 


Velocity x Slope of 


Slope of Attenuation 
" versus 

Vol. Fraction of Porosity 


-l 


(9) 


For the glass-fiber/epoxy-matrix composites used to validate our procedure the value of 
the Numerical Factor was found to be = 0.5 (using a velocity of 0.25 cm/|xsec) and for 
the graphite-fiber/epoxy-matrix composites studied the value for the Numerical Factor 
was found to be = 2 (using a velocity of 0.30 cm/|xsec). The good agreement between the 
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Correlation Plot for 



Volume Fraction of "Porosity" 


Figure 6. Correlation plot for measurements on graphite-fiber/epoxy-matrix composites. 
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experimental data and the phenomenological model suggests that this approach may be 
useful in estimating concentrations of porosity in composite laminates. 

Ill KRAMERS-KRONIG RELATIONSHIPS 
Hla. INTRODUCTION 

Previous reports from this Laboratory have focused on local approximations to the 
following Kramers-Kronig relationships 


2 7 CO'XVcu') 

X 1 (co)-X 1 (oo) = ± f 2V . d co' 

x o co - or 

( 10 ) 

2cd 7^t( co')-Ki(oo) 

/sr 2 (co) = J -= T d(Q 

X g CO 2 — CO 2 

(11) 


where K^cS) and K 2 ((£)) are the real and imaginary parts, respectively, of the dynamic 
compressibility K (co) . Tests of several approximations to these exact relations have been 
reported on previously (every progress report from the 9/15/85 to 3/14/86 report to this 
one). In the last proposal we indicated our intention to provide a derivation of the local 
approximation to Eq. (10) on which our previous investigations were based. The next 
two portions of this report. Sections IHb and IIIc, present such a derivation. They also 
demonstrate that the local approximation we have tested may be applicable in systems 
exhibiting considerable dispersion such as carbon-carbon composites. 

Rib. DERIVATION OF TWO LOCAL APPROXIMATIONS 

The nonlocal character of Equations (10) and (11) limits their usefulness. Fre- 
quency domain data extending over the range from 0 to °° is never available experimen- 
tally due to the limited bandwidth of any measurement system. Computable approxima- 
tions to Eqs. (10) and (11) are needed. The purpose of this Section is to derive a local 
approximation to the Kramers-Kronig relations. The approximation derived here is not 
the only local approximation to these equations which is possible. In fact, most of the 
common Hilbert transform pairs do not satisfy the approximation that will be derived 
here. 13 

We will present a derivation of a local approximation in a very general setting. In 
order to do this we will first rewrite the Kramers-Kronig relations in a more general nota- 
tion 

b (co) = . 

x J 0 u 2 - co 2 


( 12 ) 
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where A (co) and B (co) are the real and imaginary parts respectively of a general fre- 
quency response function 0(co). The dynamic compressibility given above K (co) is just 
one of many choices we could have made for the frequency response function. In fact, 
many researchers have taken as their starting point the acoustic index of refraction and 
then proceeded from the Kramers-Kronig relations relation between the real and ima- 
ginary parts of this quantity. This result is essentially equivalent to the approach we have 
taken as long as we assume that the attenuation coefficient is greater than zero for all fre- 
quencies but zero and =*= which seems reasonable. (This claim will be substantiated 
below, in Section IIIc.2.) The first step in the approximation scheme is make the change 
of variables 


X = m/co . 


With this change of variables the integral becomes 



A(cac)-A(co) 

x 2 -l 


dx , 


(13) 


(14) 


where the integration variable is now dimensionless. Now the first in a series of approxi- 
mations is made by writing 


1+p/G) 

B(co) = — f A ( tax )~ A (<o) dx 
n J 


l-p/G) 


1 


(15) 


where p is the distance between co and the nearest singularity of A (u ) in the complex u 
plane. If it happens that p > co then the appropriate approximation is 


„ l+p/co 

B( CO) = - J ■ 
71 0 


A (co x )-A (co) 


x 2 -l 


dx . 


(16) 


Before proceeding it is worth noting that Eq. (15) constitutes an example of the 
finite Hilbert transform. As such the exact of the approximate equation is available. 14 
However, this approach is not pursued here. Making the approximation in Eq. (15) 
involves neglecting the tails of the integral. The first term neglected is 


2 1 p / %( cojQ-Atco) ^, 
n l x 2 -l 


(17) 


which has, due to Assumption lb), an absolute value less than 


2 M 

% 


1 - p/co 


JC “1 


-dx . 


(18) 


Thus the absolute value of the first tail is less than 
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— In 
7t 


p/CO 


2 -p/co 

for p < co , and is equal to 0 for p > co. The other tail is 

2 7 A(cox)-A(co) 


^ J 


1+p/co 


x 2 -l 


dx , 


which has an absolute value that is less than 


(19) 


( 20 ) 


2 M f 1 j 

I -5 dx. (21) 

71 l+p/co* — 1 

Making the further change of variables z = x - 1 leads to 


Hence, 


2 M 7 \_ 

71 l+p/co^ 2 - 


—dx = 
1 


■J 


-dz < 


2M 


f 


-d? 


2_ 7 A (cox )-A (co) < 2Mco 

71 l+p/co JC 2 — 1 


( 22 ) 


(23) 


A review of the error bounds that have been established so far indicates that in the 
case where co « p these errors will decrease to zero. (Under some circumstances, it is 
also possible that the term that has been retained may be decreasing to zero at the same 
rate, in which case the approximation will fail.) 

The next step in the approximation scheme is to expand A (u) in an infinite series, 
using Burmann’s series expansion. 15 The expansion is given in the following theorem: 


Burmann’s Theorem: 


Let \\f(z ) be a function defined by the equation 

z - a 


V(z) = 


(24) 


<l>(z) - <Ka)’ 

where <|)(z ) is analytic in a circle centered at z=a. Then any function f(z) analytic 
in the same circle centered at a can be expanded in the form 


/w-m + ’f' w-j ai 


m = 1 


mi 


m j/n - 1 

— lf'(a)iW)} m ]+R n (25) 


da 


m 


where 


R = J_ [ dl f J <K 5 )- W 
" 27 t/J J <b(r) — <b(a ) 


1 n - 1 


<KO -<K^) ’ 


( 26 ) 
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and y is a contour in the t -plane, enclosing the points a and z and such that, if £ be 
any point inside it, the equation <KO = <K 0 has no roots on or inside the contour 
except a simple root t = £. 

Of course the choice of <]) is completely arbitrary up to the conditions imposed by the 
theorem. It is this arbitrariness which makes it possible to derive many different local 
approximations to the Kramers-Kronig relations. <|> should be chosen so that the first few 
terms of the series provide a good estimate of A (u ) over the frequency range of interest. 
This choice is extremely important to approximations based on experimental data where 
the calculation of even the first derivative is often unreliable due to noise. 

If the function A (to) is a slowly varying function of to then it would seem reason- 
able to expand A (co) in terms of a function <)> which is also varying. In this case one 
choice for <)> is 

<) >(u) = ln(u/(0) (27) 

where the argument of <J) has been chosen to be a dimensionless function of co. With this 
choice of ((), \j/ becomes 

\}/(co) = lim y(u ) = lim — “f 7 . 03 — = lim-j- = to. (28) 

u—ho u->(Dln(U/Q)) —U u-x o J_ 

U 

Inserting this into Burmann’s expansion leads to the following series representation for 

A (u) 


A(u) = A (co) + £ 


m jtn - 1 


m= 1 


{\n{u/co)} m d 
m\ dof‘ 


— [A'(co)co m ] 


(29) 


where a prime denotes differentiation with respect to co, since the expansion of A (u ) was 
carried out in the complex u -plane which contains co. Using this sum in Eq. (15) leads to 

1+p/co 

jm-1 


B (to) 


= 1 f _J_y . 
V J p/ co* 2 " 1 ^! 


{ln(jc)} m df 
m! dvf 


— [A , (co)co' n ]<ft . 


(30) 


Assuming that the series in Eq. (29) is uniformly convergent and recalling that the 
integral appearing in Eq. (30) may be thought of as a principal value integral, we see that 
it is possible to bring the summation outside of the integral to produce 


B( 0 fl = -£-1 

7C “ m ! 
m = 1 


1+p/co 

dm ^ r [A'(co)co m ] J i y) l^ 

1 — p/co X 


d (£>‘ 


- 1 


dx 


(31) 


It is still, however, appropriate to think of the integral in Eq. (10) and, hence, Eq. (30) as 
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a principal value integral. With this restriction the singularity at x = 1 is excluded from 
the integration range and thus if Burmann’s series, Eq. (29), is uniformly convergent over 
the range of integration then so is the series obtained by dividing Eq. (29) by l/(x 2 - 1).) 
The next step is to further approximate each of these integrals by moving the lower 
bound of integration back to 0 and moving the upper bound out to «*>. It is an easy matter 
to evaluate the error introduced by taking the upper bound of integration out to <*>. Error 
estimates are obtained by evaluation of 


iff* 

1+p/G) x - L 


(32) 


which goes to zero as p — » <*> for fixed co. Estimation of the error introduced by taking 
the lower bound of integration back to 0 is harder to obtain. These error estimates would 
involve evaluation of integrals like 


1— p/co 


J 


(Incor 

* 2 -l 


dx 


(33) 


all of which have a singularity at the lower bound of integration. However, if co « p 
then it will not even be necessary to evaluate these integrals since the integration range 
already extends to 0 in that case. Of course, all integrals must be taken as principal value 
integrals in order to exclude the singularity at * =0 from the range of integration. Furth- 
ermore, it is apparent that if the frequency response function is getting small around co 
then the errors introduced by moving the bounds of integration as just discussed may 
cause the approximation to fail. 

Having the range of integration as described above, we see that the expression for 
B( co) is now 


B{ CO) = - 2 ~V dm m T [A »co m ]pf {ln{ * )} -dx (34) 

k ^ ml dco 1 i x 2 -l 

The integrals in Eq. (34) may be found in many tables of integrals, at least for the first 
few values of m. 16 For the sake of completeness, however, it will be shown how to 
evaluate these integrals. The first thing to note is that these integrals vanish for even 
values of m . This may be seen as follows. Letting m = 2n , 


= pj^Ldx + ij$f¥-dx 

o x 2 — 1 J 0 ;t 2 -1 {x 2 -\ 


( 35 ) 


In the integral from 1 to °°, making the change of variables x = y gives 
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= + pfM^fLdl 

l x 2 - 1 J 1/y 2 - 1 y 

- + 


(36) 


(37) 


Now, the minus sign preceding the lny may be dropped since 2n is even. Next writing x 
in place y gives 

l l 

HnvY 2 " 

(38) 


= pf-Sf^-d* + pfi!“£-<fe = 0 

I ! -1 l 1-x 2 


Thus, the sum in Eq. (34) actually ranges over only the odd values of m . 

In order to evaluate the remaining integrals, the following contour integral will be 

used 



where the contour C is that depicted in Fig. (7) together with the subcontours which will 
be needed below. It must be stressed that the complex plane in which the contour C 
exists is not the same complex plane as was used to derive the Kramers-Kronig equa- 
tions. It might be more accurate to refer to the complex plane used in the latter case as 
the complex angular frequency plane since numbers in that domain have dimensions such 
as rad-Hz or rad-MHz whereas the complex numbers in which the contour C is placed 
are dimensionless. In order to differentiate between these two complex planes the con- 
vention will be adopted in which dimensionless complex numbers will be denoted by y 
and the set of these complex numbers referred to as the complex y plane. As a further 
convention real numbers in the complex y plane will be denoted by the symbol x . Com- 
plex numbers in the complex angular frequency plane will be denoted by u and the set of 
these complex numbers will be referred to as the complex u plane. For example, the 
expansion of A (to) carried out using Burmann’s theorem was done in the complex u 
plane and the integrals appearing as a result of the substitution of this series into the 
Kramers-Kronig equations have been rescaled, Eq. (27), so that they may be evaluated in 
the complex y plane. Note that the integral over the large semicircle C 6 vanishes as the 
radius of that semicircle goes to infinity. Additionally, the integral over the small 
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Complex y-plane 
(Dimensionless) 



Figure 7). Contour used to evaluate the integrals in Eq.(39) shown in the 
dimensionless complex y-plane. 
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semicircle C 5 of radius p centered at the origin vanishes as p approaches zero. To see 
this, parameterize the semicircle as follows. Let y = pe lQ where 0 ranges from —n to k. 
Then, 


f (lny)m dv = f OnP + fOr 

c 5 y 2 - 1 lp 2 e 2iQ -l 


pie ld dQ . 


(40) 


Using the binomial theorem permits this to be rewritten as 


m 


X 

k = 0 


m 

k 



(i Q) m ~ k ie iQ 
p 2 e 2iQ - 1 


dQ 


(41) 


Furthermore, 

— flnol^ t ) 

limp(lnp)* = lim (lnp) = lim— ^ = -)tlim-^^— ( 42 ) 

p — >0 p — >0 p ->0 1 p ->0 1 

p p 2 p 

where L’Hospital’s Rule has been used. We find, after application of this rule k times, 
that 


(-l)* -1 *!limplnp = 0 (43) 

p ->0 

So in the limit as p — >• 0 the integral over C 5 vanishes. 

We now turn to the evaluation of the contributions from the remaining contours. 
The integrals over C 4 and C 3 are evaluated using a generalization of the residue theorem. 
It must be noted, however, that the negative portion of the real axis lies on the second 
branch of the log function. That is, on the negative real axis, 

ln(y) = ln(|y|)-i7C (44) 

Bearing this in mind, we see that the contribution from C 4 is 

.jMD-ig. . tUfL ( 4S) 

(Since d(y 2 -l)/dy evaluated at y =-l is -2.) and the contribution from C 3 is 


[ln(l)] w 

2 


= 0 . 


(46) 


Now the integrals over C j and C 2 may be evaluated as principal value integrals. Collect- 
ing results from above gives 


■ p f dx +P ]SfTL dx+ ti^l. 

—00 X 1 0 X 1 


2 


( 47 ) 
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where the convention that real numbers in the complex y plane are to be denoted by x 
has been used. The binomial theorem may be used to rewrite this as 


m- 1 

o- 2 

k = 0 



m 


k 

_ 



w OO 

~ k p[^~-dx + ip [ iiHZLdx 

U 2 -\ J„x 2 -1 


(-in) 


m+1 


So 


J Onxr 

lx 2 - 1 


_ (—iitf 1 * 1 ' m 1 


■dx = 


-12 7 

fc=i L J o ^ 1 


For m = 1 the sum vanishes and 


(48) 


(49) 


jf . (Ins) 

Jx 2 -1 


dx = 


TC 


(50) 


Form =3 we get 


,f dnx) 3 
lx 2 - 1 


3 = _ (~ £ tt) 4 _ 1 2 

4 2 


dx = 


2 

*=1 


(- i n) 3 ~ k P f dx (51) 


J v 2 _ i 

0 x i 


o A 1 


_ ft 4 3 ji 7t 2 

4 4 4 

_ 3it 3 _ _7t^ 

“16 2 

Thus, from Eq. (34), the local approximation to first order is 


B( co) = - 

Jt 


Tt 2 

4 


CO 


dA (co) 


dco 


B(0» = 

2 dco 


(52) 


(53) 


To find the corresponding approximation for the other Kramers-Kronig equation, 
Eq. (11), 

oo B (m ) B (co) 


A (co) - A (0) = 


2or f u 
n l u 2 - 


co 


-du 


co 
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proceed by rewriting it in the following form 


A (co) - A (0) = - co 


B(u) fl(co) 


2co f u co 

n l u 2 - CO 2 



( 54 ) 


The expression in the parentheses is identical in form to the integral equation whose 
approximation has just been derived so long as B ( u )/u takes the same role as A (u ) took 
earlier and the expression A (co) - A (0) takes the place of B (co). Then the following is 
immediate 


A (co) — A (0) = - A (55) 

2 dco co 

An alternate local approximation may be obtained similarly by starting with 
Eq. (29) but expanding A (u ) in powers of ln(- u /co) about the point u = co, where the log- 
arithm is thought of as a complex valued function. In this case the same approximation 
scheme developed above may be used. However, the integrals evaluated in Eq. (34) 
become 


UnHOr 

6 x 2 -l 


J 


-dx 


This is equal to 

j [!"(*>->}" 

o x - 1 

For m = 1 this becomes 

j . ]MU . 1* 

0*1 0*1 0*1 

_ 7lM x)) 


_ fWii 

l V 2 1 

0 X - 1 


, K 

-dx = — . 
4 


On the other hand, the analog of Eq. (28) is 

u - co 


\|/(co) = lim 

U — 





r *s 


In 

u 

- In 

CO 



CO 

fa fal 


CO 

fa - 



= lim l/(— 1 /m ) = -co 

U — » CO 


This leads to a local approximation which is 


B( co) = f- 

K 


K 


4 


-co 


dA( co) 


d co 


_ 7t dA (co) 
2 dco 


(56) 


(57) 


(58) 


(59) 


(60) 
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and 


A (cd) - A (0) 


Tcco 2 d B (co) 
2 d co co 


This local approximation provides an accurate description of many acoustical systems 
where the attenuation coefficient is an increasing function of frequency. 17, 18 

Before going on, it is worth noting that Burmann’s expansion applied directly to 
©(co) could also be used to derive a local approximation. This could be achieved by 
expanding ©(co) in terms of some function ((> chosen to be purely imaginary or real. Tak- 
ing either the real or imaginary parts of the resulting series expansion would lead to a 
purely local relation between the real and imaginary parts of ©(co), i.e., between A (co) 
and B (co). However, in this case the right-hand- side of the relation would have terms 
containing both A (co) and B (co). The advantage obtained by proceeding as was done in 
this Section is that the relation obtained gives A (co) purely in terms of B (co) and 
vice versa. 


Bib. DISCUSSION AND PREVIEW 

There is at least one obvious way in which both local approximations may fail. 
Inspection of Burmann’s series, Eq. (26), shows that every term except the first contains 
A '(co). If A '(co) happens to be zero then the expansion technique discussed here must fail 
to give any useful information about the behavior of the function B (co). If A '(co) = 0 then 
A (co) must be a constant. If A (co) is constant globally, i.e., constant for all co, then B (co) 
must be zero also. This result follows from the global form of the Kramers-Kronig equa- 
tions. However, it is not clear a priori that if A (co) behaves locally like a constant that 
B (co) must be zero. In this sense the global and local Kramers-Kronig equations lead to 
the same predictions. 

It is also necessary to point out that there is considerable freedom in the choice of 
the expansion function <|> used in Burmann’s theorem Eq. (25). For example, the function 
ln(- u ) was used instead of ln(« ) to obtain an alternate local approximation. This led to 
the local approximation given by Eq. (60). In fact the first local approximation derived, 
Eq. (53), was tested by Booij and Thoone in viscoelastic theory and found to agree well 
with experiment. This approximation, Eq. (53), may also have application in acoustics, 
for instance, in the case where the phase velocity is decreasing with frequency but where 
the corresponding attenuation must be positive. 

It may seem surprising to approximate an integral by a derivative, especially in light 
of the fact that derivatives usually make functions rougher and integrals usually are con- 
sidered to be smoothing operators. However, the integrals analyzed in this Section were 
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singular integrals, so intuition may be misleading in this case. Indeed, the main results of 
this Section indicates that this is so. In fact, it has been shown in this Section that the 
Kramers-Kronig integrals behave more like derivatives than regular (i.e., nonsingular) 
integrals. This result has a precedent in the theory of complex variables in the form of 
Cauchy’s formula for the derivative of an analytic function expressed as a complex con- 
tour integral. Beyond that, the local character of these relations should not be surprising 
when it is recalled that the analysis began with analytic functions. The Cauchy-Riemann 
equations provide an example of a local relationship between the real and imaginary 
parts of any function at any point where that function is analytic. Results of a similar 
nature have been obtained by Widder 19 as well as Paley and Wiener 20 for the case of the 
Stieltjes transform. Examination of these results should give credibility to the idea that it 
is possible to obtain local estimates of the behavior of functions related by nonlocal 
integral equations. 

Specifically, the Stieltjes transform of <)>(r ) is 



Widder has shown that this nonlocal relation has the purely local inversion given by 


<J)(r ) = lim 


(-0 


jfc-i 


tk-i 


!(£ -2)! dt 


irfl'V (')] • 


Paley and Wiener derive an alternate local expression which is 

In 


1 m 

<K0 = lim-^-SH)" 


d 

Tit— 

dt 


(t 1/2 f(0) 


i/i ■ 

m 7 It 0 


(2n)\ 


A common feature to both of these expressions is that they call for the evaluation of 
derivatives up to infinite order. This seems to be a common feature of many local inver- 
sion formulae of nonlocal operators. 
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HIc. THE LOCAL APPROXIMATION IN ACOUSTICS 

In this Section the local approximations derived in Section Ilia, will be applied to 
acoustical systems. In Section nib.l we will derive a local approximation relating 
dispersion to attenuation in systems exhibiting small dispersion. In Section IIIc.2 we will 
generalize this relation to systems in which the dispersion is not necessarily small. Thus, 
we will fulfill one of the objectives of our previous proposal which was to generalize 
these relations to systems which exhibit considerable dispersion. In this Section we shall 
also see that the acoustic index of refraction could have been chosen as our frequency 
response function without changing any of our results. Thus our approach and that taken 
other researchers, who used the acoustic index of refraction for their frequency response 
function, are essentially equivalent. A preview of the results that we will obtain in this 
Section is presented in Table 1. 

1) The Case of Small Dispersion 

Up to now everything that has been done is perfectly general. At this point a 
specific frequency response function @(co) will be chosen, namely, the dynamic compres- 
sibility K (co) of the medium in which a traveling sound wave will be propagated. (In this 
Section of the report we will apply the local approximation, Eq. (60) to weakly dispersive 
materials. This is essentially a review of previously reported results and is included to 
provide a contrast with the material to be presented in Section IIIc.2 of the present 
report.) At this point the attenuation coefficient and phase velocity finally enter. This is 
because the frequency dependent compressibility must also satisfy the dispersion relation 
for acoustic wave propagation, 

k 2 = co 2 po/«r (co) , (61) 

where k is the wave number of the acoustic wave, and p 0 is the density of the material. 
Making the usual identification that 

‘ = ^ +i “« < 62 > 

where C( co) is the phase velocity of the acoustic wave and a(co) is the attenuation 
coefficient, and substituting this relation into Eq. (61) gives 

2 

= co 2 Po X(co). (63) 


CO 


C( co) 


+ i a(co) 


Writing K (co) as the sum of its real and imaginary parts, 

K( co) = K l (w) + iK 2 (a)), 


(64) 
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Table 1. A flow chart showing the different developments given in Section IIIc.l 
(shown in the left column) and nic.2 (shown in the right column). Both 
developments lead to the same conclusion but the path on the right is more general. 
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inserting this into Eq. (63), and equating real and imaginary parts gives 

r 1 2 


0) 

C(CQ) 

2a((o) 

C(co) 


a 2 (co) = 00^^ (co) 
■ = COpo^2(«) 


(65) 

( 66 ) 


Shortly, AT^co) will be identified with A (to) and AT 2 ( co) will be identified with B (to). 
However, before this can be done, Eqs. (65) and (66) must be decoupled in order to sim- 
plify the derivation below. This may be accomplished via the following approximation 


«(co)C(co) 

to 


« 1 


(67) 


which is true in many experimental situations. Eq. (67) states that the magnitude of the 
imaginary part of the wave vector is much less than the magnitude of the real part. 
Table 2 contains a compilation of typical values for the left-hand side of Eq. (67). The 
Table is organized in terms of increasing values of a parameter (3, the so-called slope of 
attenuation which is defined in Eq. (82). We note that Eq. (67) is very well satisfied for 
all entries in Table 2. With this assumption Eqs. (65) and (66) become 

2 


Po^t(co) = 


PoA: 2 (co) = 


1 


C (co) 
2a(co) 


cr(co) 


or 


( 68 ) 


CO 


C (co) 


Dropping the small term gives 


poAT^co) = 


C (co) 


(69) 


PoAT 2 (co) = 


2a(co) 


co 


C (co) 


These equations may now be solved for AT^co) and Af 2 (co), yielding 

1 


AT^co) = 


Po>C (co) 2 


* 2(C0) = 

p 0 coC (co) 


(70) 

(71) 


Now Eq. (60) may be applied to obtain 


«*> ■ -t^r 


(72) 



Table 2. A compilation of the ultrasonic parameters needed to compute the coupling term appearing in Eq.(162). The materials tabulated 
above cover a 20 to 1 variation in the value of the slope of attenuation, beta. The value of f was chosen to be 5 MHz for the computed 
values shown for the coupling term presented in the right-most column. 
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Substituting the expression for K ^co) in terms of C (co) gives the following expression 
for dK i(oi)/d co 


dKiiad s 1 -2 dC{ co) 

<*co p 0 c( co) 3 dco 

dC(co) __ PqC(co) 3 ^i(co) 
d co 2 d co 


Now substituting Eq. (74) into Eq. (72) gives 


dC(oi) = PqC(co) 3 2 1 
dco 2 Jt co 


■K 2 ( co) . 


The expression for K 2 ((0) in terms of a(co) may be used to obtain 


dC( co) _ 2C (co) 2 a(co) 
dco rcco 2 


(73) 

(74) 


(75) 


(76) 


from Eq. (75). This equation, Eq. (76), may be rewritten as a separable differential equa- 
tion, 


a(co) = 


rcco 2 dC (co) 
2C (co) 2 d co 


which may be solved by integration. This leads to 


2 ? aCco') ^, = ^ dC (co') = 1 

n L co' 2 ig C (co') 2 C(co 0 ) 


1 

C(co) 


(77) 


(78) 


It is worth noting that the approximation in Eq. (67) was essential to the derivation of this 
equation. It will be seen later that there is another way to derive Eq. (78) which rests 
solely on the local approximation and which is not limited by the assumption of small 
dispersion. However, at this point it is again necessary to use the fact that the material 
has small dispersion, that is. 


C(co 0 ) = C( co). 


(79) 


Then it is found that 


1 1_ C 1 (co) - C (cop) C (co) - C (co 0 ) 


C(co 0 ) C (co) C (co 0 )C (co) 

This leads finally, with C 0 = C (co 0 ), to 

2C} ® 


C (co 0 ) 


C(co)-C 0 = — J-^p-dco'. 

t r * 


71 coo CO' 


(80) 


( 81 ) 


As we have reported in previous progress reports this equation agrees well with experi- 
mental data, obtained from both homogeneous and inhomogeneous materials, for fre- 
quencies from 2 MHz to 20 MHz. However, it may hold over a much wider frequency 
range. Table 3 presents the predictions of Eq. (81) for a variety of different dependences 
of the attenuation coefficient on frequency. In certain cases further conclusions may be 
drawn from Eq. (81). For instance, experimental evidence indicates that the attenuation 
is an approximately linear function of frequency for a wide range of frequencies, i.e.. 


a(co) = -—CO 

lit 

for a wide range of materials. With this assumption, Eq. (81) predicts 

Q 2 g 

C(co) - C 0 ~ — Y-ln(co/co 0 ) . 

n 


(82) 


(83) 


HI2) Acoustical Systems where the Dispersion Is Not Necessarily Small 

As mentioned above, it is not really necessary to make the assumption that the 
dispersive effects are small. This may be accomplished by, first, showing that K (co) * 0 
in the lower half plane, which permits the square root of Eq. (63) to be taken without 
introducing any singularities in the lower half plane. This procedure produces a new sys- 
tem response function, also analytic in the lower half plane, which has real and imaginary 
parts that are functions of a(co) and C(co). The resulting Kramers-Kronig relations, 
obtained for this new system response function, are the same as those obtained by 
researchers who use the acoustic index of refraction for their frequency response func- 
tion. Before this can be demonstrated, however, a small detour into the theory of com- 
plex variables must be made. The development given here follows closely that given by 
Landau and Lifshitz. 21 It rests on two theorems from complex variables which for com- 
pleteness will be stated here. For convenience they will be called Theorem 2 and 
Theorem 3. They are: 

Theorem 2: 


Let a function f be analytic inside and on a simple closed contour C , except for at 
most a finite number of poles interior to C . Additionally, let f have no zeros on C 
and at most a finite number of zeros interior to C . Then if C is described in the 
positive sense, 



-31- 


Description 

Observed Frequency Dependence 
of Attenuation Coefficient 

Predicted Frequency Dependence 
of Phase Velocity 

Rayleigh Scattering 

<D 4 

CO 3 

Classical Viscous 

® 2 

G) 1 

Plastics, Composites 

,1 

In cg 

Soft Tissue 

CO 

Geometric Scattering 
(Hollow Pores) 

0 

(O 

(D 1 


Table 3. A compilation of the observed behavior of the frequency dependence of 
the attenuation ceofficient for several physical phenomena and the corresponding 
frequency dependence of the phase velocity as predicted by Eq.(81). 
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where N is total number of zeros off inside C and P is the total number of poles 
there. A zero of order m 0 is to be counted m 0 times, and a pole of order m p is to be 
counted m p times. 

The next theorem is known as the Argument Principle. 

Theorem 3: 


Let C be a simple closed contour described in the positive sense and let f be a 
function which is analytic inside and on C , except possibly for poles interior to C . 
Additionally, let f have no zeros on C. Then 

^-^arg f(z) = N -P (85) 

where N and P are the number of zeros and the number of poles of f , counting 
multiplicities, interior to C and A c arg/(z) is the change in the argument of f as 
we traverse the contour C once. 

If we now return to the very general setting in which the Kramers-Kronig integral 
equations were established, and, furthermore, if we suppose that @(co) is such that its 
imaginary part B (co) is never zero on the real co axis except at zero and infinity, we will 
be able show that ©(co) has no zeros in the lower half of the complex co plane. 

Our first step is to consider the integral 

1 r d_ 0(co) dco g 

lizi £ d co 0(co) - a 

where C is the contour shown in Fig. (8). According to Theorem 2 this is equal to the 
difference between the number of zeros and the number of poles of the function 
0(co) - a in the region interior to the contour C . Let a be a real number and C be the 
contour consisting of the real axis and an infinite semicircle in the lower half plane. Sup- 
pose that 0(0) = 0 O is finite. Since 0(co) is analytic in the lower half plane, so is 
©(co) - a . Therefore the integral in Eq. (86) is simply the number of zeros of the func- 
tion 0 (cd) - a , that is, the number of points at which 0 (cd) takes on the real value a . Our 
next step is to perform the conformal mapping from the co plane to the © plane via the 
transformation 

co — > © , 0 = ©(co) . (87) 

Under this transformation Eq. (86) becomes 

1 r d® 

2iti £,&-a 


( 88 ) 
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Figure 8). Drawing of the contour used in Eq.(86) and its image under the conformal 
transformation defined in Eq.(87). 
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The contour C' is the image of the contour C under the conformal transformation 
defined in Eq. (87). C ' is depicted schematically in Fig. ( 8 ). Now, the function 0(oo) 
goes to zero as co goes to » by the Riemann-Lesbegue Lemma. For this reason the entire 
infinite semicircle in the contour C is mapped onto the single point 0 = 0 on C' . The 
point co = 0 is mapped onto the point 0 = 0 O . Furthermore, 0 O is real since 
0 O = 0(0) = A (co) (since B (co) = 0 when co = 0). The right and left halves of the real co 
axis are mapped onto some very complicated (and generally self-intersecting) curves 
whose shape is not significant. The important thing about these curves is that they lie 
entirely in the upper and lower halves of the © plane respectively. This observation 
about the image of the positive co axis follows from the assumption that the imaginary 
part B (co) of 0(co) is strictly greater than zero, except at co = 0 or The corresponding 
result for the image of the negative co axis follows immediately then from the odd parity 
of the function B (co). Since B (co) is strictly greater than zero, the two images in the 
complex 0 plane nowhere intersect the real axis of the complex © plane. Because of this 
property of the curve C', it can be seen that the total change in the argument of the com- 
plex number 0 - a , as we pass once around the contour C' , is 2k if a lies between 0 and 
0 O or zero if a lies outside that range. Hence it follows from Theorem 3 that the expres- 
sion in Eq. (88) is unity for 0 < a < 0 O , and zero otherwise. 

Thus we conclude that, if a is in this range, then 0(co)-a takes the value zero once 
and only once for co in the lower half plane, but if not, then Q-a is never zero in the 
lower half plane. Hence, we deduce that on the negative imaginary co axis, where the 
function 0(co) is real, it cannot have either a maximum or a minimum since otherwise it 
would take on some values at least twice. As a result 0(co) varies monotonically on the 
negative imaginary co axis. Because 0 is analytic in the lower half plane, it can be seen 
that on the negative imaginary axis, and nowhere else, 0 takes on all real values between 
0 and 0 O once and only once. 

Let’s now return to our previous setting, and replace 0(co) by K (co). We observe 
that the imaginary part of K (co) is 

= (oporto) (89) 

cm 

Physically, it is clear that C (co) and a(co) are never less than or equal to zero for any 
positive value of co. Thus we see that K 2 (co) > 0 for all positive co. This means that 
K (co) is never zero in the lower half co plane and, as has already been shown, that it is 
analytic there. Thus the square root of K (co) must be analytic in the lower half of the co 

plane. Hence, by Eqs. (63) and (64) the physical quantity —?==• — 7 — is the real part of 

vPo C(co) 
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a Kramers-Kronig pair and that - La is its corresponding imaginary part. Inserting 

VPo o> 

these quantities into the local approximation, Eq. (60), and canceling the factor of - ■ ,L= 

VPo 

that appears on both sides of the resulting equation yields 

1 


a(oo) it d 

■ - - ~ CO 

co 2 dco C( co) 


(90) 


which leads to 


CD 

r = _i 

i o>' 2 C(coo) cm 


(91) 


This is precisely the same result that was obtained earlier in Eq. (78) by assuming that 
dispersive effects were small. However, this restriction appears nowhere in the deriva- 
tion of Eq. (91). Thus, the local approximation, Eq. (91), may hold in systems which 
exhibit considerable dispersion such as carbon-carbon composites. Table 1 reviews the 
developments given in Sections IH.c.l andIII.c.2. 
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